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Evaluation of two-loop self-energy basis integrals using differential equations 



Stephen P. Martin 
Physics Department, Northern Illinois University, DeKalb IL 60115 USA 
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I study the Feynman integrals needed to compute two-loop self-energy functions for general masses 
and external momenta. A convenient basis for these functions consists of the four integrals obtained 
at the end of Tarasov's recurrence relation algorithm. The basis functions are modified here to 
include one-loop and two-loop counterterms to render them finite; this simplifies the presentation 
of results in practical applications. I find the derivatives of these basis functions with respect 
to all squared-mass arguments, the renormalization scale, and the external momentum invariant, 
and express the results algebraically in terms of the basis. This allows all necessary two- loop self- 
energy integrals to be efficiently computed numerically using the differential equation in the external 
momentum invariant. I also use the differential equations method to derive analytic forms for various 
special cases, including a four-propagator integral with three distinct non-zero masses. 



I. INTRODUCTION 

The comparison of data with the predictions of the 
Standard Model, and candidate extensions of it, requires 
the kind of accuracy obtained from two-loop and even 
higher-order calculations. As a forward-looking exam- 
ple, if supersymmctry proves to be correct then the Large 
Hadron Collider will be able to measure the mass of the 
lightest neutral Higgs scalar boson to an accuracy of or- 
der 100 MeV, and a future linear collider will certainly do 
better 0. In contrast, even assuming perfect knowledge 
of all input parameters, the present theoretical uncer- 
tainty is probably at least 10 times larger 

The motivation for the present paper is to facili- 
tate routine calculations of self-energies, and thus pole 
masses, for particles in any field theory. A key step in 
this process is the evaluation of the necessary two-loop 
integrals. It has become clear that analytical methods 
will only work in special cases, so practical numerical 
methods are needed. In this paper, I will build on the 
many irnpqrtant advances that have been made in this 
area with the goal of streamlining both compu- 

tations and presentations of results for self-energies. 

Tarasov Q has provided a solution to the problem 
of reducing two-loop self-energy integrals to a minimal 
basis, such that any scalar integral can be represented 
as a linear combination of integrals of just four types, 
plus terms quadratic in one- loop integrals. (Other use- 
ful reduction algorithms are presented in and 0.) 
Tarasov's algorithm relies on the integration by parts 
technique Q and repeated use of recurrence relations in- 
volving integrals in different numbers of dimensions 0- 
The two-loop scalar basis integrals remaining after apply- 
ing this algorithm have the topologies shown in Figure ^ 
They are the three-propagator "sunrise" diagram S*, a di- 
agram T which is obtained from the sunrise diagram by 
differentiating with respect to one of the squared masses, 
a four-pr opag ator diagram [/, and the five-propagator 
"master" [ll| diagram M. 

Consider a generic two-loop integral Fi{s;x,y, . . .), 



which depends on the external momentum invariant 

s = -p\ (1.1) 

[using either a Euclidean or a signature ( — \ — I — h) metric] 
and propagator squared masses x,y, . . .. For special val- 
ues of the arguments, it may be possible to compute Fi 
analytically in terms of polylogarithms fsOj or Nielsen's 
generalized polvlogarithmslsij . This requires Q that 
there is no three-particle cut of the diagram for which the 
three cut masses mi, m2, m^, the invariant s for the total 
momentum flowing across the cut, and the four quantities 

s — (mi ± 7712 ± 7713)^ (1-2) 

are all non-zero. Many analytical results for various such 
special cases have been worked out There are 

also expansions [23|-[2^ in large and small values of the 
external momentum invariant, and near the thresholds 
and pseudo-thresholds ,22j-j3^]. Integral representations 
[s^-^Jl allow for systematic numerical evaluations. 

In this paper I rely instead on the differential equation 
method evaluating the basis integrals. The 
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FIG. 1; Feynman diagrams for the two-loop basis integrals. 
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idea is to take advantage of the fact that the basis inte- 
grals Fi satisfy a set of coupled first-order linear ordinary 
differential equations in s, of the form 

'i^^ = T.^vF, + a. (1.3) 

Here and Ci are ratios of polynomials in s and the 
squared masses. (If we include only genuine two-loop 
functions in the set Fi, then Ci will also include terms 
linear and quadratic in the one-loop functions, which are 
known analytically and present no problems.) The val- 
ues of the functions Fi are known analytically at s = 0. 
So one can integrate the differential equations from the 
initial conditions at s = to the desired value of s using 
well-known numerical techniques such as Runge-Kutta. 
For the integrals of the type S, T, U, this has already 
been done and explained in detail in jl^-^^. Here, I 
will extend these results to include the master integral 
M, and present results for S, T, U integrals in a different 
basis which may be more convenient for some purposes. 

In order to find the differential equations in s that the 
basis integrals satisfy, I proceed by first calculating the 
derivatives of the basis integrals with respect to their 
propagator squared-mass arguments. Using Tarasov's 
recurrence relations, these derivatives are expressed al- 
gebraically in terms of the basis functions, in the linear 
form: 

^F, = -f (1.4) 

The equations \i.'6\ in s will then follow by elementary 
dimensional analysis, using the known dependence of the 
basis functions on the renormalization scale. The deriva- 
tives of the basis functions with respect to the squared 
masses are also useful in their own right, since each 
derivative adds an extra power of the corresponding prop- 
agator in the denominator. This provides a simplified al- 
gebraic algorithm for computing integrals with arbitrary 
powers of the propagators present in the master integral 
topology. 

The rest of this paper is organized as follows. Section 
^defines the basis integrals, and gives conventions and 
notations. Section lTTll Dresents the derivatives of the basis 
integrals with respect to their squared-mass arguments. 
In section Hvl I give the differential equations in s satis- 
fied by the basis functions. The numerical integration of 



the differential equations near s = relies on expansions 
for small s, which are provided in section Section IVTI 
presents some analytic expressions for the basis functions 
in special cases that are useful both for comparison with 
the literature and for practical purposes. Section lVlIl de- 
scribes the numerical computation of the basis integrals, 
and gives two examples. 

II. CONVENTIONS AND SETUP 

The loop functions in this papers are defined by scalar 
Euclidean momentum integrals regularized by dimen- 
sional reduction to d = 4 — 2e dimensions. Let us define 
a loop factor 

C={lQ7:^)^^^{2^^lr/^\ (2.1) 

The regularization scale /i is related to the renormaliza- 
tion scale Q (in the MS scheme [s^, or the DR scheme 
[5d| for supersymmetric theories, or in the DR scheme 
|54j for softly broken supersymmetric theories) by 

q2 = 47re-V- (2-2) 

Logarithms of dimensionful quantities are always given 
in terms of 

IKX = ln(X/Q2). (2.3) 

The loop integrals are functions of a common external 
momentum invariant s as explained in the Introduction. 
(Note that the sign convention is such that for a stable 
physical particle with mass m, there is a pole at s = iv? .) 
Throughout this paper, s should be taken to have an 
infinitesimal positive imaginary part. Since all functions 
in any given equation have the same s, it will not be 
included explicitly in the list of arguments. 

The one-loop self-energy integrals [s^ are defined as: 

The two-loop integrals are defined as: 



S(x,?/,z) 
T(x,2/,z) 



[fc2 + x] [g2 + y] [[k + «7 - p)2 + z] ' 



dx 



S{x,y,z), 



1 

[fc2 + x] [{k - J3)2 + y] [g2 + z][{q + k - pf + u] 



(2.6) 
(2.7) 
(2.8) 



3 



M{x,y,z,u,v) ^ j d'^k J d'^q 



[fc2 + x] [g2 + y] [{k - p)2 + z] [{q ~ pY + u][{k ~ qY + v] 



(2.9) 



I find it convenient to introduce modified integrals in which appropriate divergent parts have been subtracted. At 
one-loop order, define the finite and e-independent integrals: 



At two loops, let 



where 



A{x) = lim [A(x) + x/ (] = x(lna; — 1), 



B[x,y) = lim [B(a;,y) - 1/e] = - / dt\n[tx + [1 - t)y - t{l - t)s 



S{x,y,z) = lim S{x,y, z) - sj^]l{x,y, z) - sj^l{x,y, z) 



(2.10) 
(2.11) 

(2.12) 



S!nl{x,y,z) = (A(x)+A(y)+A(z))/e, 
S^l{x,y,z) = {x + y + z)/2e^ + {s/2-x-y-z)/2e 

are the contributions fi'om one-loop subdivergences and froi n the remaining two-loop divergences, respectively. 



(2.13) 
(2.14) 



Also, 



T{x,y,z) = -—S{x,y,z). 



(2.15) 



Similarly, define 



U{x,y,z,u) = 

'lJ{x,y,z,u)-Uj^l{x,y)-U^l 



lim 



where 



U^l{x,y) = Bix,y)/e, 

Uj^l = -l/2e2 + l/2e 



(2.16) 



(2.17) 
(2.18) 



and, since the master integral is free of divergences, 

M{x,y,z,u,v) = lim M(x, y, z, u, w). (2.19) 

£—►0 

Thus, the bold-faced letters A,B,S,T,U represent the 
original regularized integrals that diverge as e ^ 0, while 
the ordinary letters A, 5, 5, T, C/, M are finite and inde- 
pendent of e by definition. Also, note that these inte- 
grals have various symmetries that are clear from the 
diagrams: 

• S{x, y, z) is invariant under interchange of any two 
of x,y,z. 

• T{x, y, z) is invariant under y ^ z. 

• U{x^ y, z, u) is invariant under z ^ m. 

• M{x,y, z,u,v) is invariant under the interchanges 
(x, z) <-> {y,u)^ and (x, y) <-> (z,u), and {x,y) *^ 
{u, z). 




FIG. 2: The two-loop Feynman diagram for V{x,y, z,u) 



This leads to many obvious permutations on formulas 
given below, which will not be noted explicitly. 

It is useful to define several related functions. The 
two-loop vacuum integral is 



I{x,y,z) = S{x,y,z)\s=a- 



(2.20) 



It is equal to {IQn^Y times the integral i{x,y,z) in 
and is precisely equal to the same function used in |23| ■ In 
the present paper, the analytical expression is reviewed 
in section IVII and the recurrence relation for derivatives 
in section Ivl 

The integral T(x, y, z) has a logarithmic infrared di- 
vergence as X ^ 0. This divergence must cancel from 
physical quantities, but as a book-keeping device it is 
useful to have a version of the integral r(0, x, y) with the 
infrared divergence removed: 

T(0,x,y) = hm [T(5,x,y)-f B(x,y)E5] . (2.21) 

d — >0 

Finally, for future reference we note that the topology V 
in Figure m arises quite often. When the vertical propa- 
gators are different, the result of the diagram is just the 
difference of two U functions. However, when the vertical 
propagators have the same squared mass y, it is useful to 
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define the corresponding integral 

Vix,y,z,u)^~^U{x,y,z,u). (2.22) 

In section liiil I will provide the formula expressing 
V{x,y, z,u) algebraically in terms of the other basis in- 
tegrals. 

To illustrate the usefulness of the above definitions, 
consider the most general renormalizable theory of real 
scalar fields 0i, governed by the interaction Lagrangian 

I \ijkn 

^ = -■^'mf(j)1 —(t)i(j)j(t)k ^4'i<Pj(t>k4'n- (2.23) 

Here mf, X'^'' and A*^'^" are the tree-level renormalized 
I 



masses and couplings. Then, defining the self-energy ma- 
trix function Ily (s) so that the pole masses and widths 
M, r are the solutions [H^ for complex s = — iVM of 
the eigenvalue equation 

(s-m2)5y -ny (s) =0, (2.24) 

one has: 

n».(^) = T7^n«(^) + 7T^^^\') + ■ • • ' (2-25) 



with 



n?(s) = h^^'^'^AimD-lx^'^-X^'^-Bimlml), (2.26) 

nlf (s) = -iA''^"A^"fA'='"'-A"'"-M(m2,m^,m2,m2,m2) 

+ ^A'""A-'"'"A™P^A"'"-[t/(m2,m2^,m2,m2) - U {ml ml ml ml)] / [ml - ml] 

+ 1 [X^'"'X^'''Py^P-Uimlml,mlml + i^ ^ j)] - iA^'=™"A^'='""5(m^, m,2„, ) 

+ \y'"'X^-PX''^^-PB{mlml)Bimlml + iA^-'"^™A'=™""A(m^)[A(m2) - A{mlM^l - ml] 

+ ^A'""A^-^"A™A(m2)[B(m2,m^) - ^(m^,^^)]/^ - ^1] 

+ \x^''"'X'^-PX"^-P[Iimlmlml) - I [ml, ml m^]/ [ml - ml], (2.27) 

in which the MS counterterms have been included. (Note that for degenerate masses, the function V will appear, 
as well as derivatives of the functions A,B,I.) Of course, for theories involving fermions and vectors, things are 
more complicated, but the basis functions as defined above tend to neatly organize the counterterms, at least in 
mass-independent renormalization schemes. 

In the following, a prime on a squared-mass argument of a function stands for a derivative with respect to that 
argument. This notation is particularly convenient when there are many derivatives or when some of the arguments 
are set equal after differentiation. Thus, for example, 



f{x",x,y') = lim 



f{x,z,y) 



dx^dy 

Several kinematic shorthand notations used throughout this paper are: 



(2.28) 



^xyz = x^ +y'^ + ~ 2xy ~ 2xz ~ 2yz, (2.29) 
D,^y, = s^- As^ix + y + z) + sl4{x + y + z)^ + 2A^y,] ~ s[6'ixyz + 4{x + y + z)A^y,] + Aly„ (2.30) 
A ~ s^v + s[v{v ~u — X — y — z) + {x — y){z — u)] -f (ux — yz){u + x — y — z) + v{x — z){y — u). (2.31) 



III. DERIVATIVES OF BASIS INTEGRALS WITH RESPECT TO SQUARED-MASS ARGUMENTS 

In this section, I present the results of taking derivatives of the basis integrals with respect to squared-mass 
arguments. These can be obtained straightforwardly, if tediously, from Tarasov's algorithm. The necessary recurrence 
relations have been implemented by Mertig and Scharf in the computer algebra program TARCER |49| . which was 
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used to derive or check most of the results in this section. The results below for (the equivalents of) the S and T 
functions have already been given in |45|. 

For the one-loop self-energy integral, one has: 

9 1 — — 

— B(a;, y) = [{x - y - s){B{x, y) - 2) + {x + y - s)lnx - 2y\ny] . (3.1) 

ox LAsxy 

Derivatives of the sunrise function S arc trivial, in the sense that they are already included in the basis: 

d 

— S{x,y,z) ^ -~T{x,y,z). (3.2) 

For the T function, there are two distinct derivatives. First, 
d 1 

—T{x,y,z) = — [kTxsS{x,y,z) + kTxTiT{x,y,z) + kTxT2T{y,x,z) + kTxT3T{z,x,y) + krx] , (3.3) 

ox XU sxyz 

where the coefficient functions are 

krxs = -2s^ + &s^{x + y + z) + s[2l^xy,^ -%[x'^ + y'^ + z'^)\+2{x + y + z)I^xy^ + 'i2xyz (3.4) 
krxTi = 2a;(x - s) [s^ - 2s{x + y + z) + A^.^^^ + Syz] (3.5) 

kxx = |5s^/12 -I- s'^x^x - 27/4] + s^x^x{yh\y + zhiz - 3x ~ 7y - 7z) + 51x/4 + 53{y + z)/4] 

+sx[\nx{2{z - X - y)ylmj + 2{y - x - z)z\nz + Sx^ -|- I0x{y + z) + ll(y^ + z"^) - I'iyz} - 41x^/4 
-103a;(y + z)/4 + llyz/6] -|- xhix{y\ny[{x - yf + 2z{x + y) - 3z^] + zhiz[{x - zf + 2y{x + z) - iy'^] 

-t-x(V + - 26?;z - x^) - ix^{y + z) - 5{y + z){y - zf} + ?,x^[{x + y + zf - 4(y - zf]^ 

+ {(a:^2/)} + {(x^z)}, (3.6) 

and kTxT2 is obtained from kxxTi by (x <-> y), and kxxTs is obtained from kxxTi by {x ^ z). The symmetries of the 
preceding expressions imply that 

an identity which seems somewhat remarkable since it is not immediately obvious from the symmetries of the Feynman 
diagram. When z = 0, this simplifies to: 

xS{x\ y, 0) = ySix, y", 0) = B(x, y). (3.8) 

The other derivative of the T function is given by 

d 1 

—T{x,y,z) = — [kTysS{x,y,z) + kTyTiT{x,y,z) + kTyT2T{y,x,z) + kTyT3T{z,x,y) + kTy], (3.9) 

Cfy -^sxyz 

where 

krys = -4s^ + 8s{x + y - z) + 12z'^ - 8zix + y) - A{x - yf (3.10) 
kryTi = -s^ + s^{z + 'iy- x) + s{bx^ + 6xy - iy^ - lAxz + 2yz + z^) - 'ix^{x + z) + {y - z)^ + Tx^y 

+lxz^ - bxy^ - 2xyz (3.11) 

kTyTz = -Ss'^z + %sz{x + y) + '&z'^{z-x-y) (3.12) 

kxy = \s^\^^xhiy - 2ln.T + ^] + s^[{2\ax - 3)zlnz - 3(a; + z/2)\nx\ny + {^x + &y + 2z)\nx - 20a; -|- z] 

+s[{?>x'^ -\- xy + 2xz + 3z2/2)lna;ln?/ -f Az{x ~ y — z)lna;lnz + z(10z — 4a;)lnz + 2{z'^ — 3y^ — 5x^ — Axy 
~Axz + 2yz)\nx + 41x^/2 + 2hxy/2 + llxz ~ GQz^/A] + {x^{y + z) - x^ - 5xyz + xz^ - z^/2)hia:In?/ 
+2z{y^ + z^ - ix^ + 2x{y + z) - 22;z)hia;Inz z{\^x^ - 10a;?; -I- 4a;z - 7z2)Inz -I- 2(2a;^ - ?>x^y -I- y^ 

+7x^z + 8xyz - Sy'^z - 8xz'^ + 3yz'^ - z^)hix + %{xz^ - x^ + x^y - 2>x^z - xyz + 3z^/2)| 
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+{x ^ y), 

and kTyT2 is obtained from kxyTi by {x y). For the special case of a vanishing first argument, one finds 
d 

■g^T{0,x,y) = kj^gS {0,x,y) + %^T(0, x,y) + kj^j,^T{x, Q,y) + %-r2^(y > 0, + % 



where 



krpg 



k 



TTl 

k— 

k'rjT 



l&y{y -X- s)/Aly - 



krpj, — (x y S)/I^sxy 



= S,xy{y - X - 3s)/ Al -{s + 3x + y)/Asxy 



sxy 

Syl-s"^ + sx -yx + ?/]/A^ 



xy 



2y 



A2 

sxy 



s[4:xlnx{ln.y - 1) + (4y - 8x)lny + 15a; - 7y] + {x - y)[4:xlnx{3 - Iny) + {8x + 4y)\ny 



-29x - 7y] + [s(ll/2 - 21na;) + 2ylny{lnx - 3) + (4a; - 2y)lna; - 9x + 13y] /A 



sxy 



(3.13) 



(3.14) 



(3.15) 
(3.16) 
(3.17) 
(3.18) 



(3.19) 



The derivatives of the U functions are: 



-^U{x,y,z,u) 



1 



dz 



d 



U{x,y,z,u) 



{x — y — s)U{x, y, z, u) + 2zT{z, u, x) + 2uT{u, x, z) + {3x — y + s)T{x, z, u) 
+AS{x, z, u) - 21{y, z, u) - 2{A{u) + A{x) + A{z)) + 2{x + z + u) - s/2 

\z — y — u)U (a;, y,z,u) + {u + z — y)T{z, x, u) — 2uT{u, x,z) + u + y ■ 



^yzu 



+ [{u + z- y)lnz + 2u{l - Inu) + 2y- 2z]B{x, y) 



—U {x, y, z, u) = kuuU{x, y, z, u) + kuTiT{x, z, u) + kuT2T{u, x, z) + kuTzT{z, x, u) 
oy 

+kus [Six, z,u)- {A{x) + A{z) + yl(u) + I{y, z,u))/2] + kusBix, y) + ku 
where the coefficient functions in the last expression are 



(3.20) 



(3.21) 



(3.22) 



kjju 


= [y-x- s)/Asxy + {y- z- u)/Ay^u - 1/y 




(3.23) 


kuTi 


= 2a:(s - a;)/yAs2:j; 




(3.24) 


kuT2 


= u{s-x~ y)/yAsxy +u{y + z- u)/yAy^u 




(3.25) 


kus 


= 2{s-x-y)/yAsxy 




(3.26) 


kuB 


= [{y + z — u)u\mi + {y + u— z)z\nz + {u — z)"^ — 


y^]/y^yzu 


(3.27) 


ku 


= [-s^/4 + s{z + u + 5a;/4 + y/i) -{z + u + x){x 


+ y)]/y^sxy + {U + Z- y)/Ay^u 


(3.28) 



and kuT3 is related to kuT2 by {z ^ u). Some care is needed in treating cases where the denominator Ay^u threatens 
to vanish. One finds by taking the limits that 

U{x,0,y,y') = ny,y',x)/2^T{y',y,x)/2-B{0,x)/2y (3.29) 
U{x,y,y\0) = [T{0,x,y) -T{y,0,x) - B{x,y)\ny] /2y (3.30) 
U{x,y',y,0) = -U{x,y,y',0) + i2-]^y)B{x,y'). (3.31) 

There are two types of derivatives of the master integral function M. First, 
d 

— M(a;, y, z, u, v) = kMxU\U{x, z, u, v) + kMxU2U{y, u, z, v) + kMxUsU {z, x, y, v) + kMxUiU{u, y, x, v) 



+ kMxS 



s 11 

^(a;, u, v) + S{y, z, v) + -j^B{x, z)B{y, u) - -I{x, y, v) - -I{z, u, v) 



+kMxTiT{x, u, v) + kMxT2T{y, z, v) + kMxTs.T{z, y, v) + kMxTiT{u, x, v) 
+kMxTb[T{v, X, u) + T{v, y, z)] + kMxBiB{x, z) + kMxB2B{y, u) + kMx 



(3.32) 
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where the coefficient functions are 



kMxUl 
kMxUI 
kMxU3 

kMxUi = 
kMxS 

• = 

kMxT2 = 
kMxT3 = 
kuxTi = 



kMx 



Finally, 



-u/A 

V — u 



+ s{2v — X — y — z — u) + {x — z){y — u) 



s{vx + vz + 2xz — yz — ux) + (a; — z){ux — vx + vz — yz) 



A A,,,A 



+ 



1 



Axyu A 



sv{v — X — y) — uvx + ux^ — uxy + 2vxy — vyz — xyz + y'^z 



Axyt; A 



2sv — uv + v — ux — vx + uy — vy — vz + xz — yz 



s{u — z — v) — ux + vx — uz — vz — xz + 2yz + , 



AsxzA 

xkMxs/2 - kMxUZ 

VkMxsl'2 — kMxU4 
zkMxs/"^ — kMxUl 
ukMxsf^ — kMxU2 



kMxTb = vkMxs/'^ + 



AxyvA 



s{v - X + y) + u{x + y -v) + y{x -y-2z + v) 

2 



kMxBi = -T^^u - kuxUi^'oy + {kMxTb - vkMxsl'2)\-a.v + 



-2vxy + vyz + xyz — y z 



Axyu A 



sv{x + y — v) + uvx — ux + uxy 



kMxB2 = -7-{v + z) - kMxUS^^X - kMxUl^'az + {kMxTb - vkMxs/^)^'av ~ 2kMxU4: 



+ 



A,^,A 



s{vx — UX — UZ + 2>vz + ^xz — 2yz + z"^) + {x — z)'^{u — v — z) 



-kMxs{x\5ix + ylny + zlnz + ulrni + 2v\n.v)/2 + -^{u + z + v - 2y) + 



(3.33) 
(3.34) 

(3.35) 
(3.36) 

(3.37) 

(3.38) 
(3.39) 
(3.40) 
(3.41) 

(3.42) 



(3.43) 



(3.44) 



sv{x — y — v) 



+v{ux + uy - 3xy - y"^ + 2yz) + {y - u){x - y)'^] + 



1 



2A,xzA 

—Avy — uz — 'dvz — xz — 6yz — 3z'^] + 8v'^{x — z) + x'^{3v — 3u — 5z) + z'^{u + 3u + 2x + 12y) 
—Au'^{x + z) + 4(t; + z — u)xy + 6{2y — 2u — x)vz — 4uvx — 14uxz + Auyz + 8zy^ + 3z^ 



s[4:u{u + V + y) — + 3ux — 3vx 



(3.45) 



—M{x,y,z,u,v) = kMvU\U{x,z,u,v) + kMvU2U{y,u,z,v) + kMvUd.U{z,x,y,v) + kMvU4U{u,y,x,v) 



where 



kuvUi 

kMvS 
khlvTl 
kjAvTb 
kMvBl 

kMv 



+kMvS 



s 11 

S{x, u, v) + S{y, z, v) + -^B{x, z)B{y, u) - -I{x, y, v) - -I{z, u, v) 



+kMvTlT{x, U, v) + kMvT2T{y, Z, v) + kMvT3T{z, y, v) + kMvT4T{u, X, v) 
+kMvT5[T{v, X, U) + T{v, y, z)] + kMvB\B{x, z) + kMvB2B{y, u) + kMv 



s{z — u — v) +u^ — uv + 2ux — uy + vy — uz — yz 



-21 A 

-kuvuz - x/A 

-{s + v)/ A + kMvUl + kMvU2 + kMvUS + kuvUi 

{kMvT5 + v/A)hlV - fcMi;!74lny - fcMi;!72lnU + 2s / A - 2kMvU\ - 2kMvU3 

^ [xlnx + ylny + zlnz + ulnu + 2vlnv — 2{x + y + z + u) — 5v + s/2] — kMvTb- 



(3.46) 

(3.47) 

(3.48) 
(3.49) 
(3.50) 
(3.51) 

(3.52) 
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Here kMvU2, kMvT2, kMvB2 are each respectively related to kMvUi, kMvTi, kMvBi by {x,z) ^ {y,u). Similarly, 
kMvUZi kMvT3 are each related to kMvUii kMvTi by {x,y) ^ {z,u), and kMvUi, kMvTi are related to kMvUi, kMvTi 
by {x,y) ^ {u,z). 

By repeatedly applying the identities in this section, one may obtain the results for two-loop Feynman sclf-cncrgy 
integrals with arbitrary powers of propagators in the denominator. An important example is that equations (|3.22|) 
and H3.31|l can be used to find the integral V{x, y, z, u) defined in eq. p. 22(1 and corresponding to the topology shown 
in Figure [3 



IV. DIFFERENTIAL EQUATIONS IN THE 
EXTERNAL MOMENTUM INVARIANT s 



where in each case a is summed over Q^, s, and the 
appropriate x,y, . . .. Section [llll already gave the deriva- 
tives with respect to the squared masses. The derivatives 
In this section, I present results for the derivatives of with respect to the renormalization scale are easily ob- 



the basis functions with respect to s. These are most tained from the definitions in section |II1 
easily obtained by dimensional analysis, using the facts 
that B, S, T, T, U, and M have mass dimensions 0, 2, 
0, 0, 0, and —2 respectively. Since the only dimension- 
ful quantities on which they depend are Q^, s, and the 
propagator masses, we have: 



d 



^Q^B{x,y) = 0, 
d 

x,y,z 



a 

d 

'Yo:—U{x,y,z,u) = 0, 



S2a—[sM{x,y,z,u,v)] = 0, 



(4.1) 
(4.2) 

(4.3) 
(4.4) 
(4.5) 



Q^^B(x,.) = 1, 

o d 



(4.7) 
(4.8) 



^S{x,y,z) = A{x) + A{y) + A{z) 

-.T-y-z + s/2, (4.9) 
Q'l^T{x,y,z) = -A{x)/x, (4.10) 



d 



d 



dQ- 



T{0,x,y) 



d 



Q -gQ^U{x,y,z,u) 
d_ 



Q^^k^M{x,y,z,u,v) 



1-B{x,y), 
1 + B{x,y), 
0. 



(4.11) 
(4.12) 
(4.13) 



Now, combining equations 1(3. l|l . ((4.1(1 . and 1(4.8(1 . one 
(4.6) finds 



s-^B{x,y) = 
as 



{s{x + y)- {x - yY)B{x, y) + {s - x + y)A{x) + {s + x - y)A{y) + s{x + y - s) . (4.14) 



Similarly, combining equations ((3.2(1 . 1(4.2(1 . and ((4.9(1 . one gets the result for the sunrise function 

3^5(0;, y, z) = S{x, y, z) + xT{x, y, z) + yT{y, x, z) + zT{z, x, y) - A{x) - A{y) - A{z) 
+x + y + z — s/2, 

and, from ^23, (gSJ, (g^J, and (|CTHl : 



(4.15) 



where 



s—T{x, y, z) = ctsS{x, y, z) + cttiT{x, y, z) + CTT2T{y, x, z) + ctti.T{z, x, y) + ct 
as 



(4.16) 



cts 



Cm 



D 



D 



sxyz 

1 

sxyz 



[s^ - s^{?,x + y + z) + 5(3x2 - j/2 „ ^2 _ 2xy - 2xz + IQyz) + {y + z ~ x)/\^y^ 
s^[2x + y + z\- s^l&x"^ + 3y2 + iz'^ + 3xy + 3xz + 2yz] + s[6x^ + '3y^ + 3z^ 



(4.17) 
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-5x2(y + z) - 4a;(y2 + z^) - 3yz(y + z) + 40a;yz] + x{y + z - x)A^y, - A^^^J (4.18) 
CTT2 = jr— [3s^ + s^(3z - 7a; - 5y) + s(5a;2 + - Tz^ - 6a;?; + 2a;z + Uyz) + [y + z ~ x)/\^y^] (4.19) 

^ sxyz 

Ct = jT^ — (s*[{]^x)/2-9/8] + s^[-y]^xl^y- {2y + 5x/2)]^x + yl^y + 4:3x/8 + 21y/A] 

+s^[y{x + + z)\Tixh\y — Syzlnylnz + [ixy + 9a;^/2)lna; + y(x ~ by + llz)ln?; — Tbx^ /% — 2\xy/2 
+52/^/4 — 69yz/4] + s[y{x'^ + 2xy — 3y^ ~ lOxz + 2yz + z^)lna;ln?/ + 2yz{x + 2y)lny\nz 
+ {4x'^y + xy'^ + 2y^ - 2y^z - 5xyz - 7x^/2)lnx + y{ly'^ - bx^ - 2xy + 14a;z - 22yz - %z^)\ny 
+ (57x^2 - Ix^y - Uxy^ - 37y^ + 2bxyz + 181?;2z)/4] + A,xyz[y(y ~x~ z)Ea;E?; 

+yz\ny\nz + {x^ — xy — y^ + yz)lna; + iy{x — y — z)lny — 2x^ — xy + by^ + 4j/z]| + (y z) (4.20) 

and CTT3 is obtained from ctt2 by the interchange (y <-> z). The equivalents of equations 14.15(1 and ((4.16(1 were 
found carhcr in 

For the T function, I find from equations 13. 1411 . ((4.4() . and ((4.11() . 

s^?'(0,a;,?/) = Cj^T{Q,x,y) + cifr,.{r{x,Q,y) + cifr,.^T{y,{),x) + C7pgS{{),x,y)+C7p (4.21) 

where 



[s{x + y)-{x~yY]/l\sxy (4.22) 

apj,^ = x(?,s + x + 0iy)/A,^y + 8xy[s{bx + y)-{x-yf]/Al^y (4.23) 

c^s = 2{s + x + y)/A,,y + 32sxy/Al^y (4.24) 
1 



A2 

sxy 



9s'*/8 + s''a;(lna; + 21/4) + s^a;[-32/lna;ln?/ + (lly - 5a;)lnx + 5a;/4 - 692//4] + sa;[4a;7;lna:ln?/ 



+ (7a;^ - 22xy - 9y^)lna; - 37a;^/4 + 181x?;/4] + {x - yfx[yliixlny - 3{x + y)lnx + 5x + Ay] 

+ ix^y), (4.25) 

and obtained from Cf-pj^ by x <-> y. 

The differential equation for the U function, obtained from equations 1(3.20(1 . ((3.21() . 13. 2211 . ((4.5() . and 1(4. 12() . is 

s-^U{x,y,z,u) = ([s(a; + y) ~ {x - y f]U{x,y,z,u) + x{y - a; - 3s)r(x,M, z) + [y - x ~ s)[2S'(a;, z,u) 

as Asxy V 

+uT{u, X, z) + zT(z, X, u) - I{y, z, u) + x{2 - Inx) + z(2 - Ez) + u{2 - Inu) - s/4] j . (4.26) 

The equivalent of this result was obtained earlier in • 

For the master integral M, I find from equations ((3.32() . 1(3.46(1 . ((4.6() . and 1(4.13(1 that: 

-^[sM{x,y, z,u,v)] = s [cmuiU{x, z,u,v) + CMU2U{y,u, z,v) + CMUzU{z,x,y,v) + CMUiU{u,y,x,v)] 
as 

s 11- 

+CMS S{x, u, v) + S{y, z, v) + -j^B{x, z)B{y, u) - -I{x, y, v) - -I(z, u, v) 

+cmtiT{x, u, v) + CMT2T{y, z, v) + cmt3.T{z, y, v) + cmt4T{u, x, v) 

+CMT5 [Tiv, X, u) + T(v, y,z)] + cmbiB(x, z) + CMB2B{y, u) + cm (4.27) 

where the coefficient functions are 

cmui — -r T" [s(y — X — v) + x^ + 2ux ~ vx — xy + vz — xz — yz] (4.28) 

CMS = 2t;/A - 2{cMUi + cmu2 + cmuz + cmua) (4.29) 

CMTi = a;(u + z - u)/A + xcms/'^- - 2xcmui - {x + z)cmu3 (4.30) 

CMT5 = sv/A + vcMs/^ (4.31) 

CMBi = su(lm; - 2)/A + scMt/4(2 - Iny) + sca/;72(2 - Inu) (4.32) 
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Cm = [{x — y){u — z) — (3s + x + y + z + u)v/2\ /A + cms[2u + .T + y + z + 7i — vhiv — [xlnx + ylnj/ 

+zlnz + Mlnu)/2] + [(3a; + z)cmui + (3t/ + u)cmu2 + (Sz + x)cmu3 + (3u + y)cMUi]/'2- (4.33) 

Here, the coefficient functions cmu2, cmt2, cmb2 are each respectively related to cmuii cmti, cmbi by (x, z) ^ {y, u). 
Similarly, cmus, cmt3 are each related to cmui, cmti by {x,y) ^ {z,u), and cmu4, cmta are related to cmui, cmti 
by {x,y) ^ {u,z). 



V. EXPANSIONS FOR SMALL s 



It is often useftil to have expressions for the two- loop integral functions expanded for small s. This provides the 
necessary initial data for integrating the differential equations numerically starting from s = 0. The expansions, given 
in terms of the analytically calculable vacuum function I{x, y, z), can be obtained by trying power series forms in the 
differential equations of the previous section. 

For example, for the one-loop function, one finds: 



Bix,y) 



Aiy) - Aix) 



g 

B{x,x) = — Ina; + — 



x~y 2{x-yY 
6xy{x + y)ln{y/x 

£.2 



x^ — y"^ + 2xy\n{y / x) 



6{x — y)^ 



{x - y)ix^ + y^ + lOxy) 



6x 60x^ 



(5.1) 
(5.2) 



For the two-loop functions, the most compact expressions involve derivatives of the vacuum integral. It is therefore 
useful to have a recurrence relation for taking derivatives of the vacuum function I(x, y, z): 



I{x',y,z) 



1 



^xyz 



[x-y- z)I{x, y,z) + {x-y + z)A{x)A{y)/x + (x + y ~ z)A{x)A{z)/x - 2A{y)A{z) 



+ {y + z- x)[A{x) + A{y) + A{z)] + x'^ - {y + zf 



l{x',x,0) = - (\nx- 1) /2, 
/(a:', 0,0) = -(h^a;-l)'/2-C(2). 
These follow immediately from the analysis in p^ . The function I{x, y, z) obeys 

xl{x\ y, z) + yl{x, y' , z) + zl{x, y, z) = /(x, z) - A(x) - A{y) - A{z) -^x + y + z, 

xl{x",y,z) = yl{x,y",z). 

These identities make the presentation of the following formulas quite non-unique. 
For the expansion of the sunrise integral, one finds 



—I(x", y.z) 

2 ^ ^ 8 



S{x,y,z) = I{x,y,z) + s 

^ O J LL 

^(O,^,^) = /(0,a;,a;)+s[-(h^x)/2-l/8] +sV36a;2 + 
Taking the derivative with respect to x yields 



^/(:."',y,z) + ^/(x"",y,z) 



(5.3) 
(5.4) 
(5.5) 



(5.6) 
(5.7) 



(5.8) 
(5.9) 



T(x,y,z) 



-I{x',y,z) + s 
1 



ll{x",y,z)-^I{x"',y,z) 



^Iix"',y,z)-^I{x"'\y,z)^yix""',y,z) 



T{x,0,x) = {l-liixf/2 + s/4:X + s^/72x^ + . 
The infrared-safe T function has the expansion 
1 



(5.10) 
(5.11) 



T(0,x,y) 



(x - y)' 



[{x + y)/(0, x, y) + 2A{x)A{y) - 2xA{x) - 2yA{y) + {x + yf 



11 



+ 



[4xyl{0, X, y) + {x + y){2A{x)A{y) + (x - 3y)A(x) + {y - 3x)A{y) + 4xy}] 
■ [24:xy{x + y)I{0, x, y) + \2{x + yfA{x)A{y) + {2x^ + 2()x^y - A2xy'^ - 2%y^)A{x) 



12(a;-y)6 

+ {2y^ + 2Qy'^x - A2yx^ - 28x^)A{y) - 3{x^ + y^) + 8xy{x^ + t/) + 86x'^y^] + 



T{0,x,x) = 
For the U integral, 

U{x,y,z,u) = 



It— 2 -r— 3 S r — 1 

in X — Inx — — + \ olnx + 1 



1 



y-x 



2 36x 



[I{x,z,u) - I{y,z,u)] + 



900x2 



[I5\nx - 19] 



(5.12) 
(5.13) 



2{y-x) 
x{x + y) 
'2(y-xY 



I{x" , z, u) 



[y - xf 
x{x + y) 



{I{x, z,u) — I{y, z,u)) + - — I{x',z,u) 



{y- xY 
x{x + y) 



{I{x,z,u)- I{y,z,u)) + 
(y - x)^ [y - a;)4 



lix' , z, u) 



I{x", z, u) + ^(^±4/(x"', z, u) + 



6(y - xf 



I2{y-x) 



I{x"",z,u) 



U{x,x,z,u) — 
For the master integral 
M{x,y,z,u,v) 



-I{x',z,u) + s —-I{x"',z,u) 



^-Iix"",z,u)-^I{x'"",z,u) 



(5.14) 
(5.15) 



1 s 

7 77 7 [^(a;, y, v) - I{x, u, v) - I{z, y, v) + I{z, u, v)] + — 

(X — z)[y — u) 4:(y — u)^[x — z)'' 



{[ 



Au + Az — X — y — 2v 



4x(y — m) 4:y{x~z)- 



[I{x, y, v) ~ I{x, u, v)] + {u- y)I{x, u, v) 



X — z y ~ u 

+v{x + y — v){x + y)I{x, y, v") + v{uv + vx — 2xy — 2uz + 2yz — u'^ — x'^)I{x, u, v") 
+2xy{x + y)I{x' , y' , v) — 2xu{x + u)I{x' , u' , v) + x{v — x — 3y)I{x' , y, v) 

+y{v — y — 3x)I{x, y' , v) + x{x + Ay — u — v)I{x' , u, v) + u{u + Az — x — v)I{x, u' , ii)| 



(5.16) 



M{x,y,x,u,v) = 



y-u 



[I{x',y,v) - I{x',u,v)] + s 



u + y 



u + y 



2Axyu 2{y — uY 



[I{x',y,v) - I{x',u,v)] 



ul{x' , u , v) + yl{x' , y', v)] + — ^ — '^LMI — -[u{x + u — v)I{x" , u, v) 



I2yu 

M{x,y,x,y,v) = I{x',y',v) 



2(y — u)2 ' ' ' ' 2Ayu{y — u) 

-y{x + y ~ v)I{x" , y' , v) + (v + 2u — x)I{x" ,u, v) — {v + 2y — x)I{x" , y, v)] 

[ul{x"', u, v) + yl{x", y, v)] 



2Axy 



5 + 6vl{v", X, y) + 2u(4t; -x~ y)I{v"', x,y) + v'^{v - x - y)I{v"" , x, y) 



(5.17) 



(5.18) 



In theories with massless vector bosons, special cases like M{x, y, x, 0, x) can arise, in which denominators implicit in 
the previous expressions threaten to vanish. However, those cases are easily obtained from the preceding, by noting 
that e.g. ul{x, y, u') vanishes as u — > 0, since /(x, y, u') diverges only logarithmically in that limit. 



VI. ANALYTICAL RESULTS 



As noted in the Introduction, for favorable mass and momentum configurations the basis integrals can be, and in 
many cases have been Q-H^l; computed analytically. The results for s — were given in the previous section. I 
will not consider other special values of s in this section; they do not typically arise in mass-independent (as opposed 
to on-shell) renormalization schemes. The remaining cases involve vanishing squared masses, which arise in theories 
with unbroken gauge symmetries, and as approximations to theories with large mass hierarchies. Results for these 
cases can be obtained by analytically integrating the differential equations presented in section Hvl with the initial 
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conditions of section ^ taking due care with the branch cuts. In this section, I will review results obtained in this 
manner, most of which have already been derived by dispersion relation and other methods. 
To compactify the notation, define the quantities 



They obey 



tabc — n ' ^abc — 7; • 

2a Za 



-L i^bca i^'cab ^bac f^cba -L 



These are exactly the changes of variables that occur in dilogarithm functional identities [50j , making the presentation 
of formulas below highly non-unique. To resolve branch cuts in the following consistent with the standard conventions 
for poly logarithms 50], it is crucial that s is always given an infinitesimal positive imaginary part. 
For the one-loop formulas, the well-known result is: 

B{x,y) = 2 - rsxyhix ~ tsyxlny + {All.l/s)\n{txys), (6.3) 
B{0,x) = 2-E2;+(a;/s-l)ln(l-s/a;), (6.4) 
B(0,0) = 2-h^(-s). (6.5) 

The two-loop vacuum integral is given bv |TollT5l[THTll20ll2ll45t 

I{x, y, z) = — [(x — y ~ z)ln?;lnz + (y — z — a;)lna;lnz + (z — x — y)lnxlny~\ + 2(a::lna:: + y\ny + zlnz) 



-^{x + y + z) + Al/yl 



Li2(?'xyz) + Li2(?'x2y) ~ in{rxyz) ln{rxzy) + ^ ln(j//a;) ln(z/a;) - ({2) 



(6.6) 



when X > y, z, and otherwise by the appropriate symmetry permutation of the arguments. Some special limits are 

l{0,x,y) = {x-y)[U2{y/x)-T^{x-y)\n{x/y) + {hixf/2-C{2)]+xT^x{2-]^y) + 2y]^y-5{x + y)/2, (6.7) 

l{0,x,x) = x[-]n X + Ahix - 5], (6.8) 
l{0,0,x) = a;[-(Sa;)V2-f 2lKa;-5/2-C(2)]. (6.9) 

When the masses are all very small, the two-loop basis integrals defined in this paper are 

5(0,0,0) - Hf_fh^(-s), (6.10) 

T(0,0,0) = -i [H-s)-l]\ (6.11) 

[7(0,0,0,0) = i [E(-s)-3]' + l, (6.12) 
M(0,0,0,0,0) = -6C(3)/s. (6.13) 

This should provide a useful quick comparison between other conventions and the ones used here. 
For the S and T functions with one vanishing mass and the others arbitrary, one finds T^l : 

^(O,^,^) = (y - x)[U2{txsy) + I'hirxsy)] - y(l ~ X / s)lii{txys)lnira;ys) + [{x + y + s)All.l/4s][lii{txys) - \n{rxys)] 
+ {y - a;)[Inx]^/2 - yhixhiy + {2x - s/4)In.T + (2y - s/4:)\ny + [{y'^ - a;^)/4s]ln(a::/y) 
-2a;-2y-hl3s/8, (6.14) 

T{x,0,y) = U2{txsy)+l^i2{rxsy) + ^'a{rxys)[yln{ryxs) + All.l]/s + rsyx^'a{y/x) + ^[\nx-l]^ -1. (6.15) 

Here I have deliberately chosen a presentation that does not make manifest |57j the symmetry under x ^ y. This 
makes the formulas slightly smaller, and also eases the taking of the limit y ^ 0: 

5(0,0,2;) = -xU2{s/x)-x(\^xf/2+{2x-s/2)]^x+[{x'^-.s^)/2s]\n{l-s/x) + 13s/8-[2 + Ci2)]x, (6.16) 
T(a;,0,0) = U2is/x) + (hix)'^ /2 -h^x + (l - x/s)lii{l - s/x) - 1/2 + ({2). (6.17) 
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The analytical expression for the T integral evidently cannot be obtained from those for S, T . By integrating the 
differential equation (|4.21|l . I find 

T{0:X,y) = (1 - 2isxv)Li2(ixsy) + (1 - 2r5a;y)Li2(?'xsy) + (2Ai/2/s)Li2(-a;r,,,/Ai/2) 



,1/2 



[{ln(xi,,,/Ai/2) + 21n(r,,,)}2 + (1 - Ey) ln(r,,,) + 21n(y/x) H^lil/x) 

+{5EsEy - 3E\ - 2E\ + ln(.T/y)}/2 + 2C(2)] + [{s -x- 2^lg)/s\ \n^{r.,ys) 
+ (1 - x/s) Mx/y) Mr^ys) + (1 - Ea;)[(x/s - y/s) \n{x/y) - (1 - Ey)]/2, 



(6.18) 



-Li2(s/a;) - (ln.x)V2 + Inx + (1 - x/s) ln(l - s/x){l - \nx - ln(l - s/x)} - 1/2 - C(2). (6.19) 



Useful cases for the U and V integrals that arise in unbroken gauge theories are compactly written in terms of the 
preceding integrals: 



U{x,y,y,{)) = -T(y,{),x) + (2-\ny)B{x,y) + l, 



1 TT^. 



V{x,y,y,{)) = —[T{{),x,y)-T{y,Q,x)-\nyB{x,y)\+{\ny-2)B{x,y'). 



(6.20) 
(6.21) 



The last integral was obtained using equation H3.31|l and the definition H2.22|l . Equivalent results were found in [l^ . 
Some other special limits of the U integral that can be quickly obtained using the differential equation method are: 

[/(a;, 0,0,0) = Li2(s/a;) + (l-a;/s)ln(l-s/a;)[Ea;-3 + ln(l-s/.x)] + (Ea:)V2-3Ea; + ll/2 + C(2), (6.22) 

C/(0, X, 0, 0) = (1 - x/s){U2{s/x) + [E(-s) - 2] ln(l - s/x)} - E(-s) + (Ea; - 2)72 + 7/2 + C(2), (6.23) 
[/(0,0, 0,a;) = -(1 + x/s)U2{s/x) - (^x)^ /2 - 2hix + (Ex - l)E(-s) - 2(1 - a;/s)ln(l - s/x) 

+11/2 -C(2). (6.24) 

Equivalent results were obtained in 

By integrating the differential equation H4.26|l with the first argument vanishing, I find: 

U{0,z,x,y) = (y - x)\L\2{t^sy) + l^hirxsy)]/ z - {y / s)\n{txys)\a.(rxys) + (1 - z/s)^ 

[tzxy - 1/2) Li2(l - txysTyxz) + Li2(l - TxysTyxz) " ^h{tyzx) - ^h{txzy) - ln(t2,j^3)ln(l - s/z) 
-ri{txyz,ryxs) ln(l - ) - ^{txyz^tyxs) ln(l - '^xys'^yxz ) - Vitxyz, l/txyz)[l'n-{tyzx) + In(ixzy)] 

+ {rzxy - 1/2) Li2(l - txystyxz) + Li2(l - Txystyxz) - Li2(rj^22:) - ^h{rxzy) - ln(ra;yz) ln(l - s/z) 



-il{rx 



) ln( 1 txystyxz ) 'nij'xyz ; ^yxs ) -^^(^ '^xys^yxz^ ^{^^xyz i ^ / '^xyz') ^^{j'yzx^ ~^ ^^{^'^xzy^\ 



+ {\nx + Iny - 4) ln(l - s/z)/2 + [\n[rxys) - \f / ^ + [ln(i„^.) - 1]V4 - \^^{x/y)/^ - ln(.x/y)/2} 
-l{x, y, z)/z + {^\l^y/2s)^^{txys) - Mrxys)\ - (y/z)E.xEy + {2x/z - l/2)Ex + (2y/z - l/2)Ey 
+ [(y - x)/2s\ ln(.T/y) + 5(z - .t - y)/2z + z/2s + {{y - x)/2z]E^a;, (6.25) 
where the function 

77(a, h) = ln(a6) - ln(a) - ln(6) (6.26) 

is employed to properly treat the branch cuts. As far as I know, this is the first analytical computation of a two-loop 
self-energy diagram with generic s and three distinct non-zero masses. I have checked it numerically using the method 
of the next section. 

Broadhurst has computed |llj| the master integral for the special limits needed in unbroken gauge theories: 

M(x,a;,y,y,0) - {txys) ^ {t^xs) - ^^^^^vs) - {^^^-ys) - ^{^'^1 (6-27) 
M(x,a;,0,0,0) = [i^3+(x/(x - s)) - 6C(3)] /s, (6.28) 

where 



6Li3(z) - 41n(z)Li2(z) - ln(l - z)\xi{z). 



(6.29) 
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F^{z) = 6Li3(-z)-41n(z)Li2(-z) -ln(l + z)ln2(z). (6.30) 

Another special case is 

M(a:,0,0,0,0) = [F+ix/{x-s))-F+is/is-x))-6aS)]/2s. (6.31) 

I have checked that these results are satisfied by the differential equation H4.27(l . using the other analytical results 
above. (Straightforward integration of equation (|4.27|) provides more complicated expressions, not given here, which 
are then evidently related to the above by some trilogarithm identities. The equivalence was checked numerically.) 

Some other special cases that have been computed in the literature will be omitted here for brevity. Ref. also 
found M(x,0,x,0,x), while ref. obtained the equivalent of U{x,y, 0,0) and M{x,0,y,0,0) and M(0, 0, 0, 0, x), 
and ref. [ig has M{x, 0, 0, x, 0), M{x, 0, 0, 0, x), M{x, x, x, 0), and M{x,x,x,x,0). 



VII. NUMERICAL EVALUATION BY 
DIFFERENTIAL EQUATIONS 

A method for using the differential equations in s to 
numerically compute basis integrals has been formulated 
by Caffo, Czyz, Laporta, and Remiddi in ji^-^^- We 
can now apply the same strategy to compute the values of 
all of the basis integrals, using the differential equations 
worked out in section Hvl 

Consider a master integral M{x, y, z, u, v) that occurs 
in a self-energy function. Typically, one will also need 
some or all of the basis integrals that arise from remov- 
ing one or more propagators. These can all be obtained 
simultaneously by solving the system of coupled first- 
order ordinary differential equations in the 15 dependent 
quantities 

M{x,y,z,u,v), U{x, z,u,v), U{y,u, z,v), 
U{z,x,y,v), U{u,y,x,v), T{x,u,v), T{y,z,v), 
T{z,y,v), T{u,x,v), T{v,x,u), T{v,y,z), 
S{x,u,v), S{y,z,v), B{x,z), B(y,u), (7.1) 

with X, y, z, M, V fixed and s as the independent variable. 
The relevant differential equations in addition to (|4.27|) 
are (jlT^ . l|I?T3l) . l|I?TE|l . and others obtained by 

obvious permutations. Since the B functions are known 
analytically, one need not treat them as among the de- 
pendent variables, but it is probably more economical in 
terms of computer processing time to do so. Other than 
the term involving B{x, z)B{y, u) in the differential equa- 
tion for the master integral M{x, y, z, u, v), the system of 
equations is linear. 

Standard computer numerical methods (for example, 
Runge-Kutta, or improvements thereof) are used to 
evolve the differential equations from s = to the de- 
sired s. Since the physical-sheet s is always taken to have 
an infinitesimal real imaginary part, and branch cuts lie 
along the real s axis, one should take the contour of inte- 
gration to lie in the upper-half complex plane. Reference 
^8] suggests using a rectangular contour going from to 
ih to s+ih to s+ie, where h is chosen large enough to stay 
away from singularities on the real s axis. Independence 
of the choice of h, and more generally on the choice of 
contour in the upper half-plane, provides a useful check 
on the numerical convergence. 



At the start of the contour at s ~ 0, the appearance 
of s on the left-hand sides of the differential equations 
requires that the initial data for derivatives of the basis 
functions with respect to s are provided, along with the 
initial values. (Alternatively, one can start the running at 
a point very slightly displaced from s = 0.) These are ob- 
tained from the expansions in section I find that it is 
often better to run sM(x, y, z, u, v) rather than the mas- 
ter integral itself. The method is always very fast and 
arbitrarily accurate, except sometimes when s is equal 
or extremely close to one of the thresholds where the 
denominators in the differential equations vanish. Even 
these cases can be efficiently computed without perform- 
ing special analytical expansions around the thresholds, 
as will be explained below. 

I have implemented this method in a computer pro- 
gram in order to test the method, and for use in future 
applications. When doing so, it is useful to note that 
all quantities other than s remain constant in the course 
of a Runge-Kutta routine. Therefore, although the co- 
efficients of various powers of s in the numerators and 
denominators of the coefficient functions are mildly com- 
plicated functions of x,y, z,u,v, they only need to be 
computed once. Comparison with specific numerical ex- 
amples for the master integral in ref. [40l | and the sunrise 
integrals in J7j yields agreement. Note that the first of 
these comparisons is actually a test of the equations and 
the method for all of the basis functions, not just the 
master integral M, since any error in any of the basis 
functions would feed into a discrepancy for the master 
integral. 

As an example, I consider the master integral and its 
subordinates for the case Q = 1, x = 1, y = 2, z — 3, 
M = 4, and v = 5. The result for the master inte- 
gral M(l,2,3,4,5) as a function of s is shown in Figure 
13 Although the dependence on s near the two-particle 
thresholds s = + « 7.464 and {^2 + 2)^ w 11.657 
is sharp, these points [and the three-particle thresholds 
{1 + 2 + V5)2 « 27.416 and {V2 + + V5)^ 28.970] 
did not present any numerical problems. The value of 
the master integral at s = is [/(I, 2, 5) — 7(1,4,5) — 
7(2,3,5) + 7(3,4,5)]/4 ft^ 0.28 889224. The asymptotic 
limit in which equation (|6.13|) is reasonably accurate is 
very far to the right of the end of the graph. Values for 
all of the basis integrals at s = 10 found from the simul- 
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FIG. 3: The master integral M(l, 2, 3, 4, 5), as a function of 
s. The heavier line is the real part, and the lighter is the 
imaginary part. 

taneous numerical solution to the differential equations 
are: 

M(l,2,3,4,5) = 0.71833535 + 0.39016220J 
[7(1,3,4,5) = -4.85695306- 2. 12756034^ 
[7(3,1,2,5) = -3.99263620- 1.79951451i 
[7(2,4,3,5) = -3.08641797 
[7(4,2,1,5) = -2.23235894 
S'(l,4,5) = -9.56660679 



r(4, 1,5) = 


-0.03036018 


r(5,i,4) = 


0.51591658 


r(i,4,5) = 


-3.01221172 


5(5,2,3) = 


-7.67047979 


r(5,2,3) = 


0.44677524 


r(2,3,5) = 


-1.69451693 


r(3,2,5) = 


-0.78612788 


i3(l,3) - 


0.77930384+ 1.53905980i 


5(2,4) = 


-0.05151328. 



As a second test case, consider the master integral 
M{x,0,0,x,x), which occurs in QED and QCD. This 
case does not satisfy the criterion for solvability in terms 
of generalized polylogarithms mentioned in the intro- 
duction, but a simple integral representation has been 
worked out in Following the method adopted here, 
the full system of differential equations simplifies to: 

s-^B{0,x) = — ^ [B(0,a;) 1 - s/xl, (7.3) 
as s — X 

5-^5(0,0,2:) = a;T(a;,0,0) + 5(0,0,a;)-a;h^a; 
as 

+2x-s/2, (7.4) 
s—S(x,x,x) — 3xT(x, X, x) + Six, X, x) — 3xlnx 

us , , , 

+6a;-s/2, (7.5) 



s-^T{0, 0, x) = - — \x{s - x)T{0, 0, x) + x{Ss + x)T{x, 0, 0) + 2(s + a:)5(0, 0, x) + x{s - Sx)\nx 
as (s — x)^ 

+5a;2 + 3sa;/4-9sV4], (7.6) 

s4^T(x, 0,0) = [2a;T(a;, 0,0) + 25(0,0, a;) + (s-2x)Ea; + 4x-9s/4l, (7.7) 

as s — X 

d 1 2 

s—T(x,x,x) — "7 -\2x(5s — 9x)T(x,x,x) + 2(s — 3x)S{x,x,x) — 2sx\n X 

as (s — xj{s — 9xj '- 

+ {s^ - 5sx + 18x^)]Ex - 36a;2 + 67sa;/4 - 9s^/4] , (7.8) 

s — U{x,0, X, x) — -TT \x{s — x)U{x, 0,x,x) — x{5s + 3x)T(x, x,x) — 2(s + x)S{x, x, x) 

as [s — xj 

+ (s + a;)a;(71na; — In a;- ll) + s(s + a;)/4], (7.9) 

s-^U(0,x,x,0) = [a;[7(0,a;,a;,0)-a;T(a;,0,0)-25(0,0,2:)+a;(5h^2:-h^\-7) + s/4], (7.10) 

as s — X 

4- [sM{x, 0, 0, X, x)] = - — \{s + 3x)T{x, x, x) + {s + x)T{x, 0, 0) + 25(x, x, x) + 25(0, 0, x) + sB{0, x)^ 
as (s — x) 

+ {2h{x - 4)sS(0, x) + 2a;(h^a; - 3)^ - 3s/2] . (7.11) 

I 



The value of the master integral obtained for x = 1 and as a function of s is shown in Figure 0] In this example, 
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FIG. 4: The master integral M(1,0, 0, 1, 1), as a function of 
s. The heavier line is the real part, and the lighter is the 
imaginary part. 

it turns out that there are numerical problems, but only 
extremely close to the double threshold at s = a; = 1, 
where it is known that prj 

Af (1, 0, 0, 1, 1) TT^ In 2 - 3C(3)/2 w 5.03800311. (7.12) 

In mass-independent renormalization schemes, this is not 
an issue since the tree-level mass appearing as the argu- 
ment of the function is not exactly the same as the pole 
mass where one will need to evaluate the self-energy. In 
on-shell type schemes, one could find the threshold value 
analytically, but one can also use the following general 
procedure to find threshold values with high accuracy. 
Near each threshold Sq, the loop functions have expan- 
sions of the form 

F{s) = F{so) + r [ai + 5i In r + ci In^ r] 

+r'^[a2 + b2\nr + C2ln^ r] + . . . (7.13) 

where r = 1 — s/sq- Now one can use the Runge-Kutta 
method to evaluate the loop functions at, say, several 
points s = so±n6 (for small integers n) , and then simply 
solve for the coefficients in the expansion, in particular 
F{so). In the present example, I find that choosing 6 ~ 
10~^, where there are definitely no numerical problems, 
and n — 1,2,3,4 is good enough to obtain the threshold 
values for s = x = Q = lto better than 9 significant 
digits. The results are: 

B(0, 1) = 2.00000000 

5(0,0,1) = -3.66486813 

T(1,0,0) = 2.78986813 

T(0,0, 1) = -3.78986813 

C/(0, 1,1,0) = 2.21013187 



M(1,0,0,1,1) = 5.03800311 

5(1,1,1) = -4.37500000 

T(l,l,l) = -0.50000000 
C/(l, 0,1,1) = -1.07973627. (7.14) 

Of these, the first six are checked using the analytic for- 
mulas (in3i), (Eini), iCT7ji . inii3, (Enni, and cnai, 

while the next two can now be seen "experimentally" 
to have the analytical values 5(1,1,1) = —35/8 and 
r(l, 1, 1) = -1/2 at threshold. 

To be extra-safe, a computer code can be configured 
to always trap the threshold and pseudo-threshold cases 
for evaluation in this manner. This is easy to do in an 
automated way, since the potentially dangerous points 
are always known in advance as the roots of the denomi- 
nators of the differential equations or from inspection of 
the Feynman diagrams. 



VIII. OUTLOOK 

In this paper I have studied the properties of a mini- 
mal basis of integral functions for two- loop self energies. 
These results include a complete set of formulas allowing 
for their automated numerical computation using differ- 
ential equations, following the same strategy as was put 
forward in It might be useful to review some of 

the advantages of this method: 

• The basis integrals can be computed for any values 
of all masses and s, to arbitrary accuracy. 

• All of the necessary basis integrals are obtained si- 
multaneously in a single numerical computation. 

• Branch cuts are automatically dealt with correctly 
by choosing a contour in the upper-half complex s 
plane. 

• Simple checks on the numerical accuracy follow 
from changing the choice of contour. 

The Tarasov recurrence relation algorithm [J can be 
used to reduce any two-loop self-energy to linear com- 
binations of these functions, with coefficients depending 
on the masses and couplings of the theory. Recently, I 
have used this basis and the methods of computation de- 
scribed here to obtain the leading two-loop momentum- 
dependent corrections to the neutral Higgs boson masses 
in minimal supersymmetry in a mass-independent renor- 
malization scheme. That result will appear soon. 

This work was supported by the National Science 
Foundation under Grant No. 0140129. 
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